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Abstract 

The problem of energy-efficient precoding is investigated when the terminals in the system are equipped 
with multiple antennas. Considering static and fast-fading multiple-input multiple-output (MIMO) channels, the 
energy-efficiency is defined as the transmission rate to power ratio and shown to be maximized at low transmit 
power. The most interesting case is the one of slow fading MIMO channels. For this type of channels, the 
optimal precoding scheme is generally not trivial. Furthermore, using all the available transmit power is not 
always optimal in the sense of energy-efficiency (which, in this case, corresponds to the communication-theoretic 
definition of the goodput-to-power (GPR) ratio). Finding the optimal precoding matrices is shown to be a new 
open problem and is solved in several special cases: 1. when there is only one receive antenna; 2. in the low 
or high signal-to-noise ratio regime; 3. when uniform power allocation and the regime of large numbers of 
antennas are assumed. A complete numerical analysis is provided to illustrate the derived results and stated 
conjectures. In particular, the impact of the number of antennas on the energy-efficiency is assessed and shown 
to be significant. 

Index Terms 

Energy-efficiency, MIMO systems, outage probability, power allocation, precoding. 



I. Introduction 

In many areas, like finance, economics or physics, a common way of assessing the performance of a system is 
to consider the ratio of what the system delivers to what it consumes. In communication theory, transmit power 
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and transmission rate are respectively two common measures of the cost and benefit of a transmission. Therefore, 
the ratio transmission rate (say in bit/s) to transmit power (in J/s) appears to be a natural energy-efficiency 
measure of a communication system. An important question is then: what is the maximum amount of information 
(in bits) that can be conveyed per Joule consumed? As reported in HI, one of the first papers addressing this 
issue is 13 where the author determines the capacity per unit cost for various versions of the photon counting 
channel. As shown in 0], the normalizecj^ capacity per unit cost for the well-known additive white Gaussian 
channel model Y = X + Z is maximized for Gaussian inputs and is given by limp^o ^"^^^^"^^-^ = ^ , 
where E|Xp = P and Z ~ CAA(0, cr^). Here, the main message of communication theory to engineers is that 
energy-efficiency is maximized by operating at low transmit power and therefore at low transmission rates. 
However, this answer holds for static and single input single output (SISO) channels and it is legitimate to 
ask: what is the answer for multiple-input multiple-output (MIMO) channels? In fact, as shown in this paper, 
the case of slow fading MIMO channels is especially relevant to be considered. Roughly speaking, the main 
reason for this is that, in contrast to static and fast fading channels, in slow fading channels there are outage 
events which imply the existence of an optimum tradeoff between the number of successfully transmitted bits 
or blocks (called goodput in 13] and HI) and power consumption. Intuitively, this can be explained by saying 
that increasing transmit power too much may result in a marginal increase in terms of quality or effective 
transmission rate. 

First, let us consider SISO slow fading or quasi-static channels. The most relevant works related to the 
problem under investigation essentially fall into two classes corresponding to two different approaches. The 
first approach, which is the one adopted by Verdu in m and has already been mentioned, is an information- 
theoretic approach aiming at evaluating the capacity per unit cost or the minimum energy per bit (see e.g., ||5l, 
|[6l . n, HI). In m, two different cases were investigated depending on whether the input alphabet contains or 
not a zero cost or free symbol. In this paper, only the case where the input alphabet does not contain a zero-cost 
symbol will be discussed (i.e., the silence at the transmitter side does not convey information). The second 
approach, introduced in is more pragmatic than the previous one. In Q and subsequent works lH, lITOl . 
the authors define the energy-efficiency of a SISO communication as u{p) = ^^^^^ where R is the effective 
transmission data rate in bits, rj the signal-to-noise-plus-interference ratio (SINR) and / is a benefit function 
(e.g., the success probability of the transmission) which depends on the chosen coding and modulation schemes. 
To the authors' knowledge, in all works using this approach (|[9l, H, |[TOl . ifTTI . |[T2l . |[T3l . etc.), the same 

'in (T| the capacity per unit cost is in bit/s per Joule and not in bit/J, which amounts to normalize by a quantity in Hz. 
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(pragmatic) choice is made for /: f{x) = (1 — e~"^)^, where a is a constant and N the block length in symbols. 
Interestingly, the two mentioned approaches can be linked by making an appropriate choice for /. Indeed, if 
/ is chosen to be the complementary of the outage probability, one obtains a counterpart of the capacity per 
unit cost for slow fading channels and gives an information-theoretic interpretation to the initial definition of 
|[9l . To our knowledge, the resulting performance metric has not been considered so far in the literature. This 
specific metric, which we call goodput-to-power ratio (GPR), will be considered in this paper. Moreover, we 
consider MIMO channels where the transmitter and receiver are informed of the channel distribution information 
(CDI) and channel state information (CSI) respectively. To conclude the discussion on the relevant literature, 
we note that some authors addressed the problem of energy-efficiency in MIMO communications but they did 
not consider the proposed energy-efficiency measure based on the outage probability. In this respect, the most 
relevant works seem to be |[T5l . |[T6l and lITTl . In |[T5l . the authors adopt a pragmatic approach consisting in 
choosing a certain coding-modulation scheme in order to reach a given target data rate while minimizing the 
consumed energy. In |[T6l . the authors study the tradeoff between the minimum energy -per-bit versus spectral 
efficiency for several MIMO channel models in the wide-band regime assuming a zero cost symbol in the input 
alphabet and unform power allocation over all the antennas. In lITTl . the authors consider a similar pragmatic 
approach to the one in lH, lITOl and study a multi-user MIMO channel where the transmitters are constrained 
to using beamforming power allocation strategies. 

This paper is structured as follows. In Sec. |II1 assumptions on the signal model are provided. In Sec. JIIJ the 
proposed energy-efficiency measure is defined for static and fast-fading MIMO channels. As the case of slow 
fading channels is non-trivial, it will be discussed separately in Sec. |lVl In Sec. |IVl the problem of energy- 
efficient precoding is discussed for general MIMO slow fading channels and solved for the multiple input single 
output (MISO) case, whereas in Sec. |V] asymptotic regimes (in terms of the number of antennas and SNR) are 
assumed. In Sec. |Vll simulations illustrating the derived results and stated conjectures are provided. Sec. I VII I 
provides concluding remarks and open issues. 

II. General System Model 

We consider a point-to-point communication with multiple antenna terminals. The signal at the receiver is 
modeled by: 

y(r) = H(r)x(r) + z(r), (1) 

where H is the x nt channel transfer matrix and nt (resp. n^) the number of transmit (resp. receive) 
antennas. The entries of H are i.i.d. zero-mean unit-variance complex Gaussian random variables. The vector 
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X is the nrdimensional column vector of transmitted symbols and z_ is an -dimensional complex white 
Gaussian noise distributed as A/'(0, cx^I). In this paper, the problem of allocating the transmit power between 
the available transmit antennas is considered. We will denote by Q = E[xx^] the input covariance matrix 
(called the precoding matrix), which translates the chosen power allocation (PA) policy. The corresponding 
total power constraint is 

Tr(Q) < P. (2) 

At last, the time index r will be removed for the sake of clarity. In fact, depending on the rate at which H 
varies with r, three dominant classes of channel models can be distinguished: 

1) the class of static channels; 

2) the class of fast fading channels; 

3) the class of slow fading channels. 

The matrix H is assumed to be perfectly known at the receiver (coherent communication assumption) whereas 
only the statistics of H are available at the transmitter. The first two classes of channels are considered in Sec. 
mil and the last one is treated in detail in Sec. |IV] and jV] 



III. Energy-efficient communications over static and fast fading MIMO channels 
A. Case of static channels 

Here the frequency at which the channel matrix varies is strictly zero that is, H is a constant matrix. In this 
particular context, both the transmitter and receiver are assumed to know this matrix. We are exactly in the 
same framework as lITSl . Thus, for a given precoding scheme Q, the transmitter can send reliably to the receiver 
log2 I In, + /oHQH^I bits per channel use (bpcu) with p = Then, let us define the energy-efficiency of this 
communication by: 

log2 |I„,_ +/9HQH^| 
GstaticlQ) — Tr{Q) ' ^^'^ 

The energy-efficiency Gstatic(Q) corresponds to an achievable rate per unit cost for the MIMO channel as 

defined in |[T|. Assuming that the cost of the transmitted symbol x, denoted by b[x}, is the consumed energy 

~ I (^x ; 

b[x\ = llxlP = Tv(xx^), the capacity per unit cost defined in [Tl is: Cs\ow — sup — ~ ~ . The supremum 

is taken over the p.d.f. of x such that the average transmit power is limited ]E[6[x]] < P. 
It is easy to check that: 
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Csiow = sup I , sup Iix;y) 

Q,Tr(Q)<P ^"^W) xM{xx")=Cl 

= sup _Gstatic(Q)- 

Q,Tr(Q)<P 



The second equality follows from 11181 where Telatar proved that the mutual information for the MIMO static 
channel is maximized using Gaussian random codes. In other words, finding the optimal precoding matrix which 
maximizes the energy-efficiency function corresponds to finding the capacity per unit cost of the MIMO channel 
where the cost of a symbol is the necessary power consumed to be transmitted. The question is then whether 
the strategy "transmit at low power" (and therefore at a low transmission rate) to maximize energy-efficiency, 
which is optimal for SISO channels, also applies to MIMO channels. The answer is given by the following 
proposition, which is proved in Appendix lAl 

Proposition 3.1 (Static MIMO channels): The energy-efficiency of a MIMO communication over a static 

channel, measured by Ggtatic. is maximized when Q = and this maximum is 

_ 1 Tr(HH^) 
^static -^j^ 2 ntcj^ ■ 

Therefore, we see that, for static MIMO channels, the energy-efficiency defined in Eq. ([3]) is maximized by 
transmitting at a very low power. This kind of scenario occurs for example, when deploying sensors in the ocean 
to measure a temperature field (which varies very slowly). In some applications however, the rate obtained by 
using such a scheme can be not sufficient. In this case, considering the benefit to cost ratio can turn out to be 
irrelevant, meaning that other performance metrics have to be considered (e.g., minimize the transmit power 
under a rate constraint). 

B. Case of fast fading channels 

In this section, the frequency with which the channel matrix varies is the reciprocal of the symbol duration 
(x(r) being a symbol). This means that it can be different for each channel use. Therefore, the channel varies over 
a transmitted codeword (or packet) and, more precisely, each codeword sees as many channel realizations as the 
number of symbols per codeword. Because of the corresponding self-averaging effect, the following transmission 
rate (also called EMI for ergodic mutual information) can be achieved on each transmitted codeword by using 
the precoding strategy Q : 

/2fast(Q) = Eh [log2 |ln. + pHQH^|] . (6) 
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Interestingly, iifast(Q) can be maximized w.r.t. Q by knowing only the statistics of H that is, E [HH^], 
under the standard assumption that the entries of H are complex Gaussian random variables. In practice, this 
means that only the knowledge of the path loss, power-delay profile, antenna correlation profile, etc is required 
at the transmitter to maximize the transmission rate. At the receiver however, the instantaneous knowledge of 
H is required. In this framework, let us define energy-efficiency by: 

Eh [log2|ln. + pHQH^|] 



GfastlQ) 



Tr(Q) 



(7) 



By defining g. as the i-th column of the matrix y^pHU, i G {1, . . . , nj}, U and {pi}^li an eigenvector matrix 
and the corresponding eigenvalues of Q respectively, and also by rewriting Gfast(Q) as 



Gfast(Q) = K; 



H 



log2 



In. + 



rit 



i=l 



H 



rit 

E 

i=l 



Pi 



(8) 



it is possible to apply the proof of Prop. 13.11 for each realization of the channel matrix. This leads to the 
following result. 

Proposition 3.2 (Fast fading MIMO channels): The energy-efficiency of a MIMO communication over a fast 
fading channel, measured by Gfast. is maximized when Q = and this maximum is 

1 Tr(E[HH^]) 



Gl 



fast 



In 2 



(9) 



We see that, for fast fading MIMO channels, maximizing energy-efficiency also amounts to transmitting at low 
power. Interestingly, in slow fading MIMO channels, where outage events are unavoidable, we have found that 
the answer can be different. This is precisely what is shown in the remaining of this paper. 



IV. Slow fading MIMO channels: from the general case to special cases 

A. General MIMO channels 

In this section and the remaining of this paper, the frequency with which the channel matrix varies is the 
reciprocal of the block/codeword/frame/packet/time-slot duration that is, the channel remains constant over a 
codeword and varies from block to block. As a consequence, when the channel matrix remains constant over a 
certain block duration much smaller than the channel coherence time, the averaging effect we have mentioned 
for fast fading MIMO channels does not occur here. Therefore, one has to communicate at rates smaller than 
the ergodic capacity (maximum of the EMI). The maximum EMI is therefore a rate upper bound for slow fading 
MIMO channels and only a fraction of it can be achieved (see [|271 for more information about the famous 
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diversity-multiplexing tradeoff). In fact, since the mutual information is a random variable, varying from block 
to block, it is not possible (in general) to guarantee at 100 % that it is above a certain threshold. A suited 
performance metric to study slow-fading channels |[T4l is the probability of an outage for a given transmission 
rate target R. This metric allows one to quantify the probability that the rate target R is not reached by using 
a good channel coding scheme and is defined as follows: 

Pout(Q, R) = Pr [log2 |ln. + pHQH^I < R] . (10) 

In terms of information assumptions, here again, it can be checked that only the second-order statistics of H are 
required to optimize the precoding matrix Q (and therefore the power allocation policy over its eigenvalues). 
In this framework, we propose to define the energy-efficiency as follows: 

rcn p^ - -^[^ ~ Pout(Q,fi)] .... 
m,R)- ^r(Q^ • (11) 

In other words, the energy-efficiency or goodput-to-power ratio is defined as the ratio between the expected 
throughput (see ll3l. |[20l for details) and the average consumed transmit power. The expected throughput can be 
seen as the average system throughput over many transmissions. In contrast with static and fast fading channels, 
energy-efficiency is not necessarily maximized at low transmit powers. This is what the following proposition 
indicates. 

Proposition 4.1 (Slow fading MIMO channels): The goodput-to-power ratio r(Q,i?) is maximized, in gen- 
eral, for Q 7^ 0. 

The proof of this result is given in Appendix |B] Now, a natural issue to be considered is the determination of 
the matrix (or matrices) maximizing the goodput-to-power ratio (GPR) in slow fading MIMO channels. It turns 
out that the corresponding optimization problem is not trivial. Indeed, even the outage probability minimization 
problem w.r.t. Q (which is a priori simpler) is still an open problem ifTSl . lIlTI . ll22l . This is why we only 
provide here a conjecture on the solution maximizing the GPR. 

Conjecture 4.2 (Optimal precoding matrices): There exists a power threshold Pq such that: 

• ifP<Po then Q* G arg minPout(Q, -R) ^ Q* G arg maxr(Q, iJ); 

Q Q 

• if P > Pq then r(Q, R) has a unique maximum in Q* = ^Int where p* < P. 

This conjecture has been validated for all the special cases solved in this paper. One of the main messages of 
this conjecture is that, if the available transmit power is less than a threshold, maximizing the GPR is equivalent 
to minimizing the outage probability. If it is above the threshold, uniform power allocation is optimal and using 
all the available power is generally suboptimal in terms of energy-efficiency. Concerning the optimization 
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problem associated with ([TTI ) several comments are in order. First, there is no loss of optimality by restricting 
the search for optimal precoding matrices to diagonal matrices: for any eigenvalue decomposition Q = UDU^ 
with U unitary and D = Diag(p) with p = {pi, . . . ,Pnt)< both the outage and trace are invariant w.r.t. the 
choice of U and the energy-efficiency can be written as: 

i=l 

Second, the GPR is generally not concave w.r.t. D. In Sec. IIV-B[ which is dedicated to MISO systems, a 
counter-example where it is not quasi-concave (and thus not concave) is provided. 
Uniform Power Allocation policy 

An interesting special case is the one of uniform power allocation (UPA): D = where p € [0, P] and 

rupA(p,i?) = r(^i„,,i?). 

One of the reasons for studying this case is that the famous conjecture of Telatar given in ifTSl . This conjecture 
states that, depending on the channel parameters and target rate (i.e., o"^, R), the power allocation (PA) policy 
minimizing the outage probability is to spread all the available power uniformly over a subset of ^* G {1, . . . , n^} 
antennas. If this can be proved, then it is straightforward to show that the covariance matrix D* that maximizes 
the proposed energy-efficiency function is |rDiag(e^. ), where e^. € Thus, D* has the same structure 

as the covariance matrix minimizing the outage probability except that using all the available power is not 
necessarily optimal, p* G [0, -P]. In conclusion, solving Conjecture 14.21 reduces to solving Telatar's conjecture 
and also the UPA case. 

The main difficulty in studying the outage probability or/and the energy-efficiency function is the fact that 
the probability distribution function of the mutual information is generally intractable. In the literature, the 
outage probability is often studied by assuming a UPA policy over all the antennas and also using the Gaussian 
approximation of the p.d.f. of the mutual information. This approximation is valid in the asymptotic regime of 
large number of antennas. However, simulations show that it also quite accurate for reasonable small MIMO 
systems ll23l. |[24l. 

Under the UPA policy assumption, the GPR rupA(?'5^) is conjectured to be quasi-concave w.r.t. p. Quasi- 
concavity is not only useful to study the maximum of the GPR but is also an attractive property in some scenarios 

^We denote by Se — {v€ {0, 1}"' | ^"li Vi — the set of nt dimensional vectors containing £ ones and nt — £ zeros, for all 
£€{l,...,nt}. 
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Is D* known? 


Is i (p) quasi-concave .' 


Is p* known? 


General MIMO 


Conjecture 


Conjecture 


Conjecture 


MISO 


Yes 


Yes 


Yes 


1 X 1 


Yes 


Yes 


Yes 


Large MIMO 


Conjecture 


Yes 


Yes 


Low SNR 


Yes 


Yes 


Yes 


High SNR 


Yes 


Yes 


Conjecture 



TABLE I 

Summary of proved results and open problems 



such as the distributed multiuser channels. For example, by considering MIMO multiple access channels with 
single-user decoding at the receiver, the corresponding distributed power allocation game where the transmitters' 
utility functions are their GPR is guaranteed to have a pure Nash equilibrium after Debreu-Fan-Glicksberg 
theorem |[25l . 

Before stating the conjecture describing the behavior of the energy-efficiency function when the UPA policy 
is assumed, we study the limits when p — )• and p — )• +oo. First, let us prove that lim TupaIp, ^) = 0. 

p—>-0 

Observe that lim Pout f — Int>-R) = 1 and thus the limit is not trivial to prove. The result can be proven 



p->o \ rit 

by considering the equivalent 1 + ^Tr(HH^) of the determinant 



H 



when a — ^ +oo. As the 

rit n,. 

entries of the matrix H are i.i.d. complex Gaussian random variables, the quantity Tr(HH^) = | hjj \ ^ 

j=i j=i 

is a 2nrnt Chi-square distributed random variable. Thus TupAip, R) can be approximated by: rupA(p, -R) = 

Rexp ( — - ) y — , with d = nt{2^ — l)cr^. It is easy to see that this approximate tends to zero when 

fc=0 ^ 

p — > 0. Second, note that the limit lim FupaCP) -R) = 0. This is easier to check since lim Pout ( —I, P ) = 0. 

Conjecture 4.3 (UPA and quasi-concavity of the GPR): Assume that D = ^'int- Then T\jpj^{p, R) is quasi- 
concave w.r.t. p G [Oi-f]- 

Table IIV-AI distinguishes between what has been proven in this paper and the conjectures which remain to be 
proven. 
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B. MISO channels 



In this section, the receiver is assumed to use a single antenna that is, 7ir = 1, while the transmitter can have 
an arbitrary number of antennas, nt > 1. The channel transfer matrix becomes a row vector h = {hi, hnj. 
Without loss of optimality, the precoding matrix is assumed to be diagonal and is denoted by D = Diag(p) 
with = {pi, ...,pn^). Throughout this section, the rate target R and noise level are fixed and the auxiliary 
quantity c is defined by: c = (T^(2^— 1). By exploiting the existing results on the outage probability minimization 
problem for MISO channels ll22l . the following proposition can be proved (Appendix Q. 

Proposition 4.4 (Optimum precoding matrices for MISO channels): For a// ^ E {1, — 1}, let be the 

|2 



unique solution of the equation (in x) Pr 



< X 



i=l 



Pr 



< X 



1=1 



where X,: are i.i.d. 



+00, c„,j = 0. Let z/„j be the 



zero-mean Gaussian random variables with unit variance. By convention cq 

nt — l i 

unique solution of the equation (in y) — — = 0. Then the optimum precoding matrices have the 

following form: 

p 



D* 



f Diag(e,) if P G 



(13) 



where c = a^{2^ — 1) and G Si. 

Similarly to the optimal precoding scheme for the outage probability minimization, the solution maximizing 
the GPR consists in allocating the available transmit power uniformly between only a subset i < nt antennas. 
As i.i.d entries are assumed for H, the choice of these antennas does not matter. What matters is the number 
of antennas selected (denoted by i), which depends on the available transmit power P: the higher the transmit 
power, the higher the number of used antennas. The difference between the outage probability minimization 
and GPR maximization problems appears when the transmit power is greater than the threshold — — . In 
this regime, saturating the power constraint is suboptimal for the GPR optimization. The corresponding sub- 
optimality becomes more and more severe as the noise level is low; simulations (Sec. IVll) will help us to 
quantify this gap. 

Unless otherwise specified, we will assume from now on that UFA is used at the transmitter. This assumption 
is, in particular, useful to study the regime where the available transmit power is sufficiently high (as conjectured 
in Proposition 14.11 ). Under this assumption, our goal is to prove that the GPR is quasi-concave w.r.t. p G [0, P] 
with D = and determine the (unique) solution p* which maximizes the GPR. Note that the quasi-concavity 
property w.r.t. p is not always available for MISO systems (and thus is not always available for general MIMO 
channels). In Appendix a counter-example proving that in the case where = 1 and rit = 2 (two input 
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single output channel, TISO) the energy -efficiency T^iso (Diag(^),i?) is not quasi-concave w.r.t. p = {pi,P2) 
is provided. 

Proposition 4.5 (UPA and quasi-concavity (MISO channels)): Assume the UPA, Q = ;^Int. then r{p,R) 
is quasi-concave w.r.t. p G [O, P] and has a unique maximum point in p* = min | -i)"tg- ^ p| ^j^gj-g j^^^ 
the solution (w.r.t. y) of: 



y 



nt-1 A 

K-l)! ^i! 



(14) 



k=l 



Proof: Since the entries of h are complex Gaussian random variables, the sum ^^|/ifcp is a 2nt— Chi- 
square distributed random variable, which implies that: 

i?{l-Pr[log2 (l + ^h^'h) <R]] 



■-MISO 



{p,R) 




(15) 



nt — 1 



R X e p 



j=0 



1 



with d = cnt = (2^ - 1) n^cr''. The second order derivative of the goodput R 



e p 



i=0 



W.r.t. p is 



R 



J^A.^-d/p^d-{nt + l)p) 



. Clearly, the goodput is a sigmoidal function and has a unique inflection 
point in pQ = ^^qr^-- Therefore, the function r^^^^(p, i?) is quasi-concave |[26l and has a unique maximum in 
p* = min|^,p| where f„j is the root of the first order derivative of T^^^'^{p, R) that is, the solution of 
d. ■ 
The SIMO case {nt = 1, > 2) follows directly since |I + pphh^\ = 1 + pph^h- 

To conclude this section, we consider the most simple case of MISO channels namely the SISO case (n^ = 1, 
rir = 1). We have readily that: 

rSisO(p,fi) = £li. (16) 
p 

To the authors' knowledge, in all the works using the energy-efficiency definition of U for SISO channels, 
the only choice of energy-efficiency function made is based on the empirical approximation of the block error 



rate which is 



■, M being the block length and x the operating SINR. Interestingly, the function given 



by (fT6l) exhibits another possible choice. It can be checked that the function e ? is sigmoidal and therefore 
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pSiSO jg quasi-concave w.r.t. p 1261 . The first order derivative of r^^^° is 

9p P'^ 



(17) 



The GPR is therefore maximized in a unique point which p* = c = cr^(2^ — 1). To make the bridge between 
this solution and the one derived in lH for the power control problem over multiple access channels, the optimal 
power level can be rewritten as: 



where E|/ip = 1 in our case. In HI, instantaneous CSI knowledge at the transmitters is assumed while here 
only the statistics are assumed to be known at the transmitter. Therefore, the power control interpretation of 
(fTSl ) in a wireless scenario is that the power is adapted to the path loss (slow power control) and not to fast 
fading (fast power control). 



In this section, we first consider the GPR for the case where the size of the MIMO system is finite assuming 
the low/high SNR operating regime. Then, we consider the UFA policy and prove that Conjecture 14.31 claiming 
that FupaIPi R) is quasi-concave w.r.t. p (which has been proven for MISO, SIMO, and SISO channels) is also 
valid in the asymptotic regimes where either at least one dimension of the system {nt, n^) is large but the SNR 
is finite. Here again, the theory of large random matrices is successfully applied since it allows one to prove 
some results which are not available yet in the finite case (see e.g., |[T9l . |[28l for other successful examples). 

A. Extreme SNR regimes 

Here, all the channel parameters {iit, ^^r> and P in particular) are fixed. The low (resp. high) SNR regime 
is defined by a'^ — )• +oo (resp. o"^ — )• 0). In both cases, we will consider the GPR and the optimal power 
allocation problem. 

1) Low SNR regime: Let us consider the general power allocation problem where D = Diag(p) with 
P = (Pii • • • :Pnt)- In ||22|, the authors extended the results obtained in the low and high SNR regimes for 
the MISO channel to the MIMO case. In the low SNR regime, the authors of |[22l proved that the outage 
probability Pout(Diag(p), i?) is a Schur-concave (see |[29l for details) function w.r.t. p. This implies directly 
that beamforming power allocation policy maximizes the outage probability. These results can be used (see 
Appendix E]) to prove the following proposition: 




(18) 



V. Slow fading MIMO channels in asymptotic regimes 
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Proposition 5.1 (Low SNR regime): When cr^ — )■ +00, the energy-efficiency function r(Diag(p), R) is Schur- 
concave w.r.t. p and maximized by a beamforming power allocation policy D* = PDiag(ei). 

2 } High SNR regime: Now, let us consider the high SNR regime. It turns out that the UPA policy maximizes 
the energy-efficiency function. In this case also, the proof of the following proposition is based on the results 
in |[22i (see Appendix Ell. 

Proposition 5.2 (High SNR regime): When cr^ — t- 0, the energy-efficiency function T(Diag{p), R) is Schur- 
convex w.r.t. p and maximized by an uniform power allocation policy D* = ^Int with p* G (0, P]. Furthermore, 
the limit when p — t- such that ^ ^ ^ is T ^^I^^, R^ — t- +00 which implies that p* — t- 0. 

In other words, in the high SNR regime, the optimal structure of the covariance matrix is obtained by 
uniformly spreading the power over all the antennas, D* = ^I^t the same structure which minimizes the 
outage probability in this case. Nevertheless, in contrast to the outage probability optimization problem, in 
order to be energy-efficient it is not optimal to use all the available power P but to transmit with zero power. 



B. Large MIMO channels 

The results we have obtained can be summarized in the following proposition. 

Proposition 5.3 (Quasi-concavity for large MIMO systems): If the system operates in one of the following 
asymptotic regimes: 

(a) nt < +00 and rir — )• +00; 

(b) rit — )■ +00 and rir < +00; 

(c) rit — )• +00, rir — ^ +00 with lim — = f3 < +00, 

rii— >+oo,«e{t,r} rit 

then rupA(j') R) quasi-concave w.r.t. p G [0, P]. 

Proof: Here we prove each of the three statements made above and provide comments on each of them 
at the same time. 

Regime (a): rit < +00 and — )• 00. The idea of the proof is to consider a large system equivalent of the 
function TupaIPi^)- This equivalent is denoted by T'^-pp^{p, R) and is based on the Gaussian approximation 



of the mutual information logg I + ^HH^ (see e.g., |[30l ). The goal is to prove that the numerator of 
rup^(p, R) is a sigmoidal function w.r.t. p which implies that V^-pxiP^ ^) ^ quasi-concave function |[26l . In 
the considered asymptotic regime, we know from |[30l that: 



log2 



I + ^HH^ 

nt 



M [nt log2 ( 1 + —pp 



nr 



■nt 



nr 



log2(e) 



(19) 
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A large system equivalent of the numerator of rupA(p, R), which is denoted by Na{p, R), follows: 



Naip,R) = RQ 



R-ntlog^il + ^pp] 



(20) 



where Q{x) - ^ 
of Na{p, R) w.r.t. p 



1 r+oo 



dt. Denote the argument of Q in (l20l) by aa- The second order derivative 



d^Na{p,R) 
dp^ 



[aa{p){a'a{p)f - oi'aip)] exp 



aa{pf 



Therefore Na {p, R) has a unique inflection point 



Pa 



rirP 



. 3/2 



1 



(21) 



(22) 



Clearly, for each equivalent of FupaCp, R), the numerator has a unique inflection point and is sigmoidal, which 
concludes the proof. In fact, in the considered asymptotic regime we have a stronger result since lim pa = 0, 
which implies that Na{p,R) is concave and therefore r^pj^{p, R) is maximized in p* = as in the case of 
static MIMO channels. This translates the well-known channel hardening effect |[30l . However, in contrast to 
the static case, the energy-efficiency becomes infinite here since T\jpa{p, R) — > | with — )• 0. 

Regime (b): nt — )■ +oo and Ur < +oo. To prove the corresponding result the same reasoning as in (a) is 
applied. From |f30l we know that: 

' PP 



l0g2 



I + ^HH^ 




nt 





■M nrlog2{l + pp), 



nt 1+ pp 



(23) 



(24) 



A large system equivalent of the numerator of rupA(p, R) is Ni,{p, R) = RQ (a(,(p)) with 

= \ — log2(e) ''' ^ ^^ [R-Hr log2(l + pp)]- 

V nr PP 

The numerator function Ni,{p, R) can be checked to have a unique inflection point given by: 

Pb = a" (2^ - 1) (25) 

and is sigmoidal, which concludes the proof. We see that the inflection point does not vanish this time (with 
nt here) and therefore the function A^;,(p, R) is quasi-concave but not concave in general. From ||26l . we know 
that the optimal solution represents the point where the tangent that passes through the origin intersects the 
S-shaped function RQ {ah{p))- As nt grows large, the function Q {ah{p)) becomes a Heavyside step function 
since Vp < pb, lim„j^_|.oo Q (afe(p)) = and Vp > pb, lim^j^+oo Q (afe(p)) = 1- This means that the optimal 
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power pI that maximizes the energy-efficiency approaches pb as nt grows large, p^ — )■ ct^ ^2"'- — 1^ The 
optimal energy-efficiency tends to ^''^p^^^) _^ i_ — when nt — > +oo. 



Regime (c): nt — )• +oo, nr — )• oo. Here we always apply the same reasoning but exploit the results derived 
in EH. From ED, we have that: 



log2 



I + ^HH^ 

nt 



M(ntni,aj) (26) 



where /// = /31og2(l + pp{l - 7)) - 7 + log2(l + PpW - 7)), o"/ = - log2 (l - l 



7 = l(l + /3+i- /(l + /3 + - 4/5 ) . It can be checked that {a'^{p)f adp) - a"{p) = has a unique 



2 \^ ' ' ' pp y y ' t-^ ' pp 
solution where a,(p) = We obtain a'^{p) = and 

a'^{p) = ^^^^^^ ^^^^^ equation {a',{p)f a^p) - 

a'^{p) = 0, there are terms in n^, n^, nt and constant terms w.r.t. nt. When nt becomes sufficiently large the 
first order terms can be neglected, which implies that the solution is given by pi{p) = 0. It can be shown 
that /x/(0) = and that m is an increasing function w.r.t. p which implies that the unique solution is p^ = 0. 
Similarly to regime (a) we obtain the trivial solution j»* = 0. 



VI. Numerical results 

In this section, we present several simulations that illustrate our analytical results and verify the two conjec- 
tures stated. Since closed-form expressions of the outage probability are not available in general, Monte Carlo 
simulations will be implemented. The exception is the MISO channel for which the optimal energy-efficiency 
can be computed numerically (as we have seen in Sec. IIV-BI ) without the need of Monte Carlo simulations. 

UFA, the quasi-concavity property and the large MIMO channels. 

Let us consider the case of UPA. In Fig. [T] we plot the GPR Pupa (p, R) as a function of the transmit 
power p G [0, P] W for an MIMO channel where nj. = nt = n with n G {1, 2, 4, 8} and p = 10 dB, R = I 
bpcu, P = 1 W. First, note that the energy-efficiency for UPA is a quasi-concave function w.r.t. p, illustrating 
Conjecture 14.31 Second, we observe that the optimal power p* maximizing the energy-efficiency function is 
decreasing and approaching zero as the number of antennas increases and also that Pupa {p* ,R) is increasing 
with n. In Fig. |2j this dependence of the optimal energy-efficiency and the number of antennas n is depicted 
explicitly for the same scenario. These observations are in accordance with the asymptotic analysis in subsection 
IV-B] for Regime (c). 



November 23, 2010 



DRAFT 



16 



Similar simulation results were obtained for the case where nt is fixed and Ur is increasing, thus illustrating 
the asymptotic analysis in subsection IV-BI for Regime (a). 

In Fig. [3l we plot the energy-efficiency Fupa {p, R) as a function of the transmit power p G [0, W for 
MIMO channel such that = 2, nj G {1, 2, 4, 8} and p = 10 dB, = 1 bpcu, P = 1 W. The difference w.r.t. 
the previous case, is that the optimal power p* does not go to zero when nt increases. This figure illustrates the 

R 

results obtained for Regime (b) in section IV-BI where the optimal power allocation pi = 0.0414 W 

and the optimal energy-efficiency Frrp^ — )• ^ = 12,07 bit/Joule when nt — +oo. 

2(2 ".--1) 

UFA and the finite MISO channel 

In Fig. 151 we illustrate Proposition 14.41 for nt = 4. We trace the cases where the transmitter uses an optimal 
UFA over only a subset of £ 6 {1, 2, 3, 4} antennas for p = 10 dB, = 3 bpcu. We observe that: i) if P < ^ 



c c 

Cl C2 



then 



c c 

C2 C3 



then the beamforming PA is the generally optimal structure with D* = P Diag(e^); ii) if P G 
using UFA over three antennas is the generally optimal structure with D* = P/2 Diag(e2); iii) if P G 
then using UFA over three antennas is generally optimal with D* = P/3 Diag(e3); iv) if P > ^ then the UFA 
over all the antennas is optimal with D* = ^ min -P| l4- The saturated regime illustrates the fact that it 
is not always optimal to use all the available power after a certain threshold. 
UFA and the finite MIMO channel 

Fig. |5] represents the success probability, 1 — Pout (D, in function of the power constraint P for nt = nr = 
2, R = 1 bpcu, p = 3 dB. Since the optimal FA that maximizes the success probability is unknown (unlike the 
MISO case) we use Monte-Carlo simulations and exhaustive search to compare the optimal FA with the UFA 
and the beamforming FA. We observe that the result is in accordance with Telatar's conjecture. There exists a 
threshold 5 = 0.16 W such that if P < 6, the beamforming FA is optimal and otherwise the UFA is optimal. Of 
course, using all the available power is always optimal when maximizing the success probability. The objective 
is to check whether Conjecture 14.21 is verified in this particular case. To this purpose. Fig. |6] represents the 
energy-efficiency function for the same scenario. We observe that for the exact threshold 6 = 0.16 W, we 
obtain that if P < 5 the beamforming FA using all the available power is optimal. If P > 6 the UFA is optimal. 
Here, similarly to the MISO case, we observe a saturated regime which means that after a certain point it is 
not optimal w.r.t. energy-efficiency to use up all the available transmit power. In conclusion, our conjecture has 
been verified in this simulation. 

Note that for the beamforming FA case we have explicit relations for both the outage probability and the 
energy-efficiency (it is easy to check that the MIMO with beamforming FA reduces to the SIMO case) and 
thus Monte-Carlo simulations have not been used. 
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VII. Conclusion 

In this paper, we propose a definition of energy-efficiency metric which is the extension of the work in [H to 
static MIMO channels. Furthermore, our definition bridges the gap between the notion of capacity per unit cost 
Ifn and the empirical approach of IH in the case of slow fading channels. In static and fast fading channels, the 
energy-efficiency is maximized at low transmit power and the corresponding rates are also small. On the the 
other hand, the case of slow fading channel is not trivial and exhibits several open problems. It is conjectured 
that solving the (still open) problem of outage minimization is sufficient to solve the problem of determining 
energy-efficient precoding schemes. This conjecture is validated by several special cases such as the MISO case 
and asymptotic cases. Many open problems are introduced by the proposed performance metric, here we just 
mention some of them: 

• First of all, the conjecture of the optimal precoding schemes for general MIMO channels needs to be 
proven. 

• The quasi-concavity of the goodput-to-power ratio when uniform power allocation is assumed remains to 
be proven in the finite setting. 

• A more general channel model should be considered. We have considered i.i.d. channel matrices but 
considering non zero-mean matrices with arbitrary correlation profiles appears to be a challenging problem 
for the goodput-to-power ratio. 

• The connection between the proposed metric and the diversity-multiplexing tradeoff at high SNR has not 
been explored. 

• Only single-user channels have been considered. Clearly, multi-user MIMO channels such as multiple 
access or interference channels should be considered. 

• The case of distributed multi-user channels become more and more important for applications (unlicensed 
bands, decentralized cellular networks, etc.). Only one result is mentioned in this paper: the existence of a 
pure Nash equilibrium in distributed MIMO multiple access channels assuming uniform power allocation 
transmit policy. 



November 23, 2010 



DRAFT 



18 



n 


= 8 


-+- n 


= 4 


n 


= 2 


— n 


= 1 



^60 




Transmit power p [W] 



Fig. 1. Energy-efficiency (GPR) vs. transmit power p e [0,1] W for MIMO cliannels wliere rir = nt = n G {1,2,4,8}, UPA D = ^Int, 
p = 10 dB, R = 1 bpcu. Observe that the energy-efficiency is a quasi-concave function w.rt. p. The optimal point p* is decreasing and FypA (P* i ^) 
is increasing vifith n. 



120r 



■3 




'l 23456789 10 



Number of antennas n 



Fig. 2. Energy-efficiency vs. the number of antennas n for MIMO rir = nt = n G {1, 2, 4, 8}, UPA, D = ^Int, p = 10 dB, B. = 1 bpcu and 
P = 1 W. Observe that FypA (p* , R) is increasing with n. 
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p = 58 mW 
,Jp',R)=12,76 biVJoule 
= 64 mW 
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Transmit power p [W] 
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Fig. 3. Energy-efficiency vs. transmit power p e [0, 1] W for MIMO Ur = 2, m € {1, 2, 4, 8}, UFA D = ^Int, p = 10 dB, R = 1 bpcu. Observe 
tiiat the energy-efficiency is a quasi-concave ffinction w.rt. p. The optimal point p* is not decreasing with n but almost constant. 
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c/c =0.6363 W 
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Power constraint P [W] 
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X 1 = 1 
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Fig. 4. Optimal energy-efficiency vs. constraint power for MISO nt = 4, rir = 1, UFA over a subset of ^ £ {1,2,3,4} antennas, p = 10 dB. 
R = 3 bpcu. We illustrate the results of Proposition 14.41 If P < ^ is low enough, the beamforming PA with full power is optimal. If P > ^ is high 
enough, the UFA is optimal but not with full power necessarily ^p* = min{^, P}^ which explains the saturated regime. 
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Power constraint P [W] 
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Fig. 6. Optimal energy-efficiency vs. power constraint P, comparison between beamforming PA, UPA and General PA for MIMO nt = rir = 2, 
R = 1 bpcu, p = 3 dB. We observe that our Conjecture I4.2l is validated. For the exact same 8 = 0.16 W, we have that for P < 5 the beamforming 
PA structure optimal and above it, UPA structure is optimal. 
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Appendix A 
Proof of Proposition 13. II 

As Q is a positive semi-definite Hermitian matrix, it can always be spectrally decomposed as Q = UDU^ 
where D = Diag(pi, . . . ,p„ J is a diagonal matrix representing a given PA policy and U a unitary matrix. Our 
goal is to prove that, for every U, Gstatk is maximized when D = Diag(0, 0, 0). To this end we rewrite 

^static as 



log2 



GstaticCU Diag(pi, . . . ,p„J U^) 



nt 



1=1 



i=l 



(27) 



where g. represents the i column of the x nt matrix G = -^/pHU and proceed by induction on rij > 1. 
First, we introduce an auxiliary quantity (whose role will be made clear a little further) 



nt 



.4=1 



log2 



i=l 



(28) 



and prove by induction that it is negative that is, V(pi, . . . ,p„J G R"', E^'^*\pi, . . . ,Pnt) < 0. 



For nt = 1, we have E^^\pi) = Tr (In^+pig^g^) ^Pig^g^ 
derivative of E^^\pi) w.r.t. pi is: 

dEW 



log2 



dpi 



-Pi[gf{In.+Pig^gf)-\^?<0 



. The first order 



(29) 



and thus E^^^pi) < E^^^O) = 0. 

Now, we assume that E^^'*~^\p) < and want to prove that E^'^*\p,pnt) < 0, where p = (pi, . . . ,p„j„i) 
It turns out that: 



dPnt 



Y.P3 



-1 



H 



— nt 



1=1 



< 0, 



(30) 



and therefore . . . ,p„,-i,PnJ < E^ipi, . . . ,p„,_i,0) = . . . ,pn,^i) < 0. 

As a second step of the proof, we want to prove by induction on > 1 that 

arg max G^"^^^ {p, Pnt ) = 0- 
For = 1 we have G).^i^^^^^ {pi ) = — =i=3_ = which reaches its maximum in pi = 0. 



(31) 
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Now, we assume that argmax G^J^l^^p{p) = and want to prove that arg max = 0. 



Let k = arg min Tr 

ie{l,...,nt} 

w.r.t. one obtains that: 



-1 



. By calculating the first order derivative of G, 



in.) 
static 



static 



dpk 



2 ' 



(32) 



with 



\i=l 
-log2 



and thus < § - ^ 



(33) 



2=1 

and Pfc = for all pi, . . . • • • ,Pnf We obtain that 



. . . ,pfc-l,0,Pfc+l, • • • 

= F'''^*~^\pi:- ■ ■ ,Pk~i,Pk+i:- ■ ■ ,Pnt), which is maximized when (pi, . . . . . . ,p„J = by 

assumption. We therefore have that Q* = UOU^ = is the solution that maximizes the function Gstatic(Q)- At 
last, to find the maximum reached by Gstatic one just needs to consider the the equivalent of the log2 |ln^ + pHQH^I 
around Q = 



log2 +pHQH 



H\ P 



Tr(HH^) 

nt 



(34) 



and takes Q = ^I^^ with q ^ 0. 



Appendix B 
Proof of Proposition 14. II 

The proof has two parts. First, we start by proving that if the optimal solution is different than the uniform 
spatial power allocation P* / ^Int with p G [O, P] then the solution is not trivial P* / 0. We proceed by 
reductio ad absurdum. We assume that the optimal solution is trivial P* = 0. This means that when fixing 
{p2, ■ ■ ■ ,Pnt) = (0, • • • , 0) the optimal pi € [0, P] that maximizes the energy-efficiency function is = 0. The 
energy-efficiency function becomes: 

,1-Pr [log2(l + m||/iif) <i?] 



T{Dmg{pi,0,...,0),R) = R- 



Pi 



(35) 
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where hi represents the first column of the channel matrix H. Knowing that the elements in are i.i.d. 



n,. 



hij ~ CA/'(0, 1) for allj £{!,... ,nr} we have that ~ expon(l). The random variable = | hij \ ^ 

i=i 

is the sum of n,. i.i.d. exponential random variables of parameter A = 1 and thus follows an 2nr chi-square 

rir — l h 

distribution (or an Erlang distribution) whose c.d.f. is known and given by qix) = 1 — exp(— x) — . We 
can exphcitly calculate the outage probability and obtain the energy-efficiency function: 

r(Diag(pi,0,...,0),ii) =iiexp (--) J] (36) 



where c = ^" ^ > 0. It is easy to check that lim T{pi,R) = 0, lim r{pi,R) = 0. By evaluating the first 

P Pi— >0 Pi— s>oo 

derivative w.r.t. pi, it is easy to check that the maximum is achieved for p\ = > where is the unique 
positive solution of the following equation (in y): 

Considering the power constraint the optimal transmission power is pi = min{ ^J'"^ , P}, which contradicts 
the hypothesis and thus if the optimal solution is different than the uniform spatial power allocation then the 
solution is not trivial P* / 0. 

Appendix C 
Proof Proposition 14.41 

Let = (pi,...,p„J be the vector of powers allocated to the different antennas i G {!,..., nf} and thus 
D = Diag(p). Define the two sets: C{x) = <p>0, ^^Pi < x> and A{x) = < p > 0, ^^Pi = x >. Using 
these notations, they key observation to be made is the following: 



£eC(P) peC(P) ^ 



sup r^isO(D,i?) R sup 



Tit 

Pi 

<i> R su, sup "-^.T°(D./i) (38) 

x6[o,P]peA(x) X 

= it sup — 

xe[o,P] ^ 



where P^^^ = Pr 



nt 



log [l + p^pi\h,\^ <R 



i=l 



: (a) translates the definition of the GPR; (b) follows from the 



property sup{AU B} = sup{sup{j4}, supjS}} for two sets A and B, applied to our context; in (c) the function 
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9[z 



z e 



) = {giiz), 



Ci-i ' Ce 



is a piecewise continuous function where gi{z) = 1 — Pr 



i=l 



for 



and £ E {1, . . . , nj}. The function g{z) corresponds to the solution of the minimization problem 



of the outage probability II221 . 

Now, we study the function g^. By calculating the first order derivative of ^.gi (^) w.r.t. x we obtain: 



d 



1 



j=0 



(39) 



Thus the function j^g (^) is increasing for x G (0, x^) and decreasing on x G (rE^,oo). The maximum point is 

jj) = where 

1 1 



reached in X£ = ^ where y£ is the unique positive solution of the equation 



We have that (f>{0) = -1 < and 



(40) 



j=0 



i=0 

i - 1 



(41) 



i=0 

> 0. 



This implies that ye < £ and thus X£ > c. Since c„ i > 1 we also have X£ > — - — for all £ G {1, . . . , — 1}. 



Therefore, all the functions (f ) are increasing on the intervals ( 



Moreover, on the interval 



oo ) , they are increasing on 



and decreasing on [x£,oo). Proposition 14.41 follows directly. 



Appendix D 
Counter-example, TISO 

Consider the particular case where nt = 2 and rir = I. From Proposition |4.4[ it follows that for a 
power constraint P < the beamforming power allocation policy maximizes the energy-efficiency and 
rTisO(Diag(p,o),i?) = rTisO(Diag(o,p),ij) > pTiso j^piag (^|, |) . The function rTisO(Diag(pi,P2), i?) 
with (pi,P2) £ ^2 — {(^1,^2) G I Pi +P2 < -P} denotes the energy-efficiency function. We want to prove 
that r^^^^(Diag(pi,p2)5 ^) is not quasi-concave w.r.t. {pi,P2) G '^2- This amounts to finding a level 7 > 
such that the corresponding upper-level setU-y = {{pi,p2) & ^2 \ r^^^^(Diag(pi,p2), ^) > 7} is not a convex 
set (see |[32l for a detailed analysis on quasi-concave functions). Consider an arbitrary < q < min|p, ^| 
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such that rTisO(Diag(g o),i?) = rTisO(Diag(o, g), ii) < T^iso (Djag |) It turns out that all upper- 
level sets U^^ with 7^ = r^^^*-'(Diag((?, 0), i?) are not convex sets. This follows directly from the fact that 
(g, 0), (0, q) G U^^ but (f , f ) ^ U^^ since pTiso (oj^g (|, f ) , i?) < 7,. 

Appendix E 
Extreme SNR cases, GPR 

In II22I . the authors proved that in the low SNR regime the outage probability Pom\.{p,R) is Schur-concave 
w.r.t. p. This means that for any vectors p, q such that p )^ q then Poat{p,R) < Po\it{(l,R)- The operator >- 



denotes the majorization operator which will be briefly described (see II29II for details). For any two vectors 

mm rit rit 

p,q & M"', p majorizes q (denoted hy p >~ q) if ^Pfc > ''^Qk, for all m G {1, . . . , — 1} and ^Pfc = ^9fc- 

k=l k=l k=l k=l 

This operator induces only a partial ordering. The Schur-convexity and -< operator can be defined in an analogous 
way. Also, an important observation to be made is that the beamforming vector majorizes any other vector, 
whereas the uniform vector is majorized by any other vector (provided the sum of all elements of the vectors 
is equal). Otherwise stated, xe^ >- p >- ^1 for any vector p such that ^^Pi = x and 1 = (1, 1, . . . , 1) and 

i=l 

Ci G Si. 

It is straightforward to see that if i^out(Diag(p), i?) is Schur-concave w.r.t. p then 1 — i^out(Diag(p), i?) is 
Schur-convex w.r.t. p. Since the majorization operator implies the sum of all elements of the ordered vectors 
to be identical, r{Diag{p), R) = ^ ^^^^ j^g Schur-convex w.r.t. p and thus is maximized by 



i=l 



a beamforming vector. Using the same notations as in Appendix O we obtain: 

sup_ r(Diag(p),i2) = sup_ - sup [1 - Pout(Diag(p), i?)] 

p6C(P) x-e[0,P] ^ peA{a:) 



(a) 



sup -[1- Pr[log(l + xphfh^) < R], 



x-e[o,P] ^ 

sup — 1 — Pr 

xe[o,P] ^ 

5n„ 



1 



li, |2 



(42) 



sup — 

a;e[0,P] ^ 



where (a) follows by considering beamforming power allocation policy on the first transmit antenna (with no 
generality loss) and replacing p = xe^ with e^^ = (1, 0, . . . , 0) and hi denoting the first column of the channel 
matrix; in (c) we make use the definition in Appendix O for the function ^Qn^ ( ) which has a unique 
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optimal point in min|^,p|, with the unique solution of ^nXv) — 0- Since cr^ 
thus the optimal power allocation is p* = Pe^. 
Similarly, for the high SNR case we have: 



then c — )• +oo and 



sup r(Diag(p),i?) 

peC(P) 



sup - sup [1 - Pout(Diag(p),i?)] 

a:6[0,7^ ^ peA(x) 



1 

sup — 

X6[0,P] ^ 



out 



Diag -(1, 

nt 



We have used the results in 11221 . where the UPA was proven to minimize the outage probability. 
Let us now consider the limit of the energy-efficiency function when p — 0, cr^ — such that ^ 



^ a positive finite constant. We obtain that 1 — 



out 



Pr 



In +^HH 

"■r rit 



H 



(43) 



with 



> which implies 



directly that F 



-In,, R 



+ 00. 



References 

[I] S. Verdu, "On channel capacity per unit cost", IEEE Trans, on Inf. Theory, vol. 36, no. 5, pp. 1019-1030, Sep. 1990. 

[2] J. R. Pierce, "Optical channels: Practical limits with photon counting", IEEE Trans. Conimun., vol. 26, pp. 1819-1821, Dec. 1978. 
[3] M. Katz, and S. Shamai, "Transmitting to colocated users in wireless ad hoc and sensor networks", IEEE Trans, on Inf. Theory, 

vol. 51, no. 10, pp. 3540-3562, Oct. 2005. 
[4] D. J. Goodman, and N. Mandayam, "Power Control for Wireless Data", IEEE Personal Communications, vol. 7, pp. 48-54, Apr. 

2000. 

[5] A. El Gamal, M. Mohseni, and S. Zahedi, "Bounds on capacity and minimum energy-per-bit for AWGN relay channels", IEEE 

Trans, on Inf. Theory, vol. 52, no. 4, pp. 1545-1561, Apr. 2006. 
[6] X. Cai, Y. Yao, and G. Giannakis, "Achievable rates in low-power relay links over fading channels", IEEE Trans, on Communications, 

vol. 53, no.l, pp. 184-194, Jan. 2005. 
[7] Y. Yao, X. Cai, and G. Giannakis, "On energy-efficiency and optimum resource allocation in wireless relay transmissions", IEEE 

Trans. Wireless Communications, vol. 4, no. 6, pp. 2917-2927, Nov. 2005. 
[8] A. Jain, S. R. Kulkami, and S. Verdu, "Minimum energy per bit for Gaussian broadcast channels with common message and 

cooperating receivers", Proc. Forty-Seventh Annual Allerton Conference on Communication, Control, and Computing, Monticello, 

USA, Sep. 2009. 

[9] V. Shah, N. B. Mandayam ,and D. J. Goodman, "Power control for wireless data based on utility and pricing", IEEE Proc. of the 
9th Intl. Symp. Personal, Indoor, Mobile Radio Communications (PIMRC), Boston, MA, pp. 1427-1432, Sep. 1998. 

[10] C. U. Saraydar, N. B. Mandayam, and D. J. Goodman, "Efficient power control via pricing in wireless data networks", IEEE 
Trans, on Communications, vol. 50, No. 2, pp. 291-303, Feb. 2002. 

[II] F. Meshkati, H. V. Poor, S. C. Schwartz, and N. B. Mandayam, "An energy-efficient approach to power control and receiver design 
in wireless data networks", IEEE Trans, on Comm., vol. 53, no. 11, pp. 1885-1894 , Nov. 2005. 

[12] S. Buzzi and H. V. Poor, "Joint receiver and transmitter optimization for energy-efficient CDMA communications", /. Set. Areas 
in Comm., vol. 26, no. 3, pp. 459-472, Apr. 2008. 



November 23, 2010 



DRAFT 



27 



[13] S. Lasaulce, Y. Hayel, R. El Azouzi, and M. Debbah, "Introducing hierarchy in energy games", IEEE Trans, on Wireless 

Communications, vol. 8, no. 7, pp. 3833-3843, Jul. 2009. 
[14] L. H. Ozarow, S. Shamai, and A. D. Wyner, "Information theoretic conisderations for cellular mobile radio", IEEE Trans, on 

Vehicular Technology, vol. 43, no. 2, pp. 359-378, May 1994. 
[15] S. Cui, A. J. Goldsmith, and A. Bahai, "Energy-efficiency of MIMO and cooperative MIMO techniques in sensor networks", IEEE 

Journal on Selected Areas in Communications, vol. 22, no. 6, pp. 1089-1098, Aug. 2004. 
[16] S. Verdu, "Spectral efficiency in the wideband regime", IEEE Trans, on Inf. Theory, vol. 48, no. 6, pp. 1319-1343, Jun. 2002. 
[17] S. Buzzi, H. V. Poor, and D. Saturnino, "Energy-efficient resource allocation in multiuser MIMO systems: A game-theoretic 

framework", Proc. of 16th European Signal Processing Conference (Eusipco), Lauzanne, Switzerland, Aug. 2008. 
[18] E. Telatar, "Capacity of multi-antenna gaussian channels", European Transactions on Telecommunications, vol. 10, no. 6, pp. 

585-596, Nov/Dec. 1999. 

[19] L. Zheng, and D. N. C. Tse, "Diversity and multiplexing: A fundamental tradeoff in multiple-antenna channels", IEEE Trans, on 

Inf. Theory, vol. 49, no. 5, pp. 1073-1096, May 2003. 
[20] S. Shamai, and I. Bettesh, "Outages, expected rates and delays in multiple-users fading channels", Proc. Conf. Information Sciences 

and Systems (CISS), Princeton, NJ, USA, Mar. 2000. 
[21] M. Katz, and S. Shamai, "On the outage probability of a multiple-input single-output communication link", IEEE Trans, on 

Wireless Comm., vol. 6, no. 11, pp. 4120-4128, Nov 2007. 
[22] E. A. Jorswieck, and H. Boche, "Outage probability in multiple antenna systems", European Transactions on Telecommunications, 

vol. 18, pp. 217-233, 2006. 

[23] Z. Wang, and G. B. Giannakis, "Outage mutual information of space-time MIMO channels", IEEE Trans, on Inform. Theory, vol. 

50, no. 4, pp. 657-662, Apr. 2004. 
[24] A. L. Moustakas, S. H. Simon, and A. M. Sengupta, "MIMO capacity through correlated channels in the presence of correlated 

interferers and noise: A (not so) large N analysis", IEEE Trans, on Inform. Theory, vol. 49, no. 10, pp. 2545-2561, Oct. 2003. 
[25] D. Fudenberg, and J. Tirole, "Game Theory", MIT Press, 1991. 

[26] V. Rodriguez, "An Analytical Foundation for Ressource Management in Wireless Communication", IEEE Proc. of Globecom, San 

Francisco, CA, USA, pp. 898-902, , Dec. 2003. 
[27] L. Zheng, and D.N. Tse, "Optimal Diversity-Multiplexing Tradeoff in Multiple Antenna Channels", Proc. Allerton Conf. Comm., 

Control, Computing, Monticello, pp. 835-844, Oct. 2001. 
[28] J. Dumont, W. Hachem, S. Lasaulce, P. Loubaton, and J. Najim, "On the capacity achieving covariance matrix of Rician MIMO 

channels: an asymptotic approach", IEEE Trans, on Inform. Theory, Vol. 56, No. 3, pp. 1048-1069, Mar. 2010. 
[29] A. W. Marshall, and I. Olkin, "Inequalities: Theory of majorization and its applications". New York: Academic Press, 1979. 
[30] B. M. Hochwald, T. L. Marzetta, and V. Tarokh, "Multiple-antenna channel hardening and its implications for rate feedback and 

scheduling", IEEE Trans, on Inform. Theory, vol. 50, no. 9, pp. 1893-1909, Sep. 2004. 
[31] M. Debbah, and R. R. Miiller, "MIMO channel modeling and the principle of maximum entropy", IEEE Trans, on Inform. Theory, 

Vol. 51, No. 5, pp. 1667-1690, May 2005. 
[32] S. Boyd, and L. Vandenberghe, "Convex Optimization", Cambridge University Press, 2004. 



November 23, 2010 



DRAFT 



